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We extend the idea of a discrete-time quantum walk on a graph by placing a qubit on each vertex, 
and allowing the walker to interact with the qubit at its current position. We show that allowing 
for a controlled-Z interaction at each time step between the walker and the vertex qubits generates 
multipartite entanglement between the qubits. We demonstrate that for particular coin operators the 
system generates a highly entangled cluster state, of use in measurement based quantum computation. 
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I. INTRODUCTION 

Quantum walks m, the quantum analogue of classical 
random walks, are conceptually simple protocols that 
display a rich variety of quantum behaviours. Both 
discrete-time |2| and continuous-time variants ISj 
can perform universal quantum computation fH [Q. 
Furthermore, they display Anderson-like localization ||6], 
and indeed this fact can be used to implement a 
quantum algorithm for search on a spatial lattice |I7|, 
offering quadratic speedup over the best known classical 
algorithms. 

Quantum walk based schemes for the generation and 
distribution of entanglement are not new |8l. However, 
the current literature mostly focuses on generating 
entanglement between the walker's external degrees of 
freedom and the coin driving its dynamics or, in two- 
dimensional quantum walks, the two motional states of 
a given walker 0. 

In this paper we present a conceptually new way 
of using quantum walks to generate entangled states, 
which we dub the interacting quantum walk. In this 
paradigm, the walker interacts with a collection of qubits, 
located at the sites of a general lattice, by means of a 
controlled-Z interaction. Using an extensive computer 
simulation of the d 5 mamics of both the lattice of qubits 
and the walker, we show that this scheme is effective in 
generating multipartite entanglement between the qubits 
in the walker itself. Moreover, interesting and useful 
forms of entangled states involving the lattice qubits 
alone can be produced by arranging suitable interactions 
with the hovering walker. 

Experimentally, considerable success has been found 
in linear-optics implementations of quantum random 
walks based on integrated photonics circuits. Recent 
advances in this context have allowed for the 


* JL's current address: Department of Computer Science, University 
College London, London, WCIE 6BT, UK 


management of disparate degrees of freedom of photonic 
information carriers, which would allow for the 
implementation of controlled interactions at the basis 
of our proposal. Our scheme, thus, holds the potential 
to drive the experimental research in quantum random 
walks well into the domain of entanglement generation 
and distribution. 

The remainder of this paper is organised as follows: in 
Sec. 1^ we introduce the formalism for the description 
of quantum walks in discrete time, and present the 
protocol for multipartite entanglement distribution that 
we have designed. Sec. Ill is dedicated to the presentation 
of the main results of our computer simulations, 
and the characterisation of the quality of the states 
produced through the proposed method. Finally, in 
Sec. 113 we draw our conclusions and motivate further 
investigations. 


II. BACKGROUND 

A. Discrete-time quantum walks 

The state space of a quantum walk on a graph G = 
{V,E) is a bipartite Hilbert space Hw = Ti-p ® 'He 
consisting of a position space 

"Hp := span({|z; e U)}) = (1) 

and a coin space 


•He := span({|d e [A])}) ^ C^, (2) 

where A = max^gy [deg(v)] is the maximum degree of 
the graph. The quantum walk evolves by successive 
alternating applications of a coin operator C : He 
He, followed by a shift operator S : Hw —^ Hw, 
which moves the walker to another vertex of the graph, 
dependent on the coin state. 

In this paper we consider quantum walks on path 
graphs L„ and cycle graphs Cn- These topologies utilize 
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FIG. 1. Illustration of the scheme of principle of the situation 
we consider. A quantum walker moves across the sites of a 
given lattice according to the 'coin tossing' operation defined 
in the main text. The walker interacts with the 'internal' degree 
of freedom of the particle occupying a given lattice site. By 
hovering over the lattice, the walker brings about information 
on the state of the lattice particles, and entangles them in a 
multipartite fashion. In panel (a) [(b)] we show the process on 
a path graph Loo [cyclic graphs Goo]. 


the shift operators 


:= \n - 1) (n - l|p 0 |0) (1|<. + |0) (0|p ® |1) (01^ 



(3) 


n—1 

.=0 ^ (4) 

+ |(f + 1) mod n) (f|p (g) |1) (l|c), 


respectively. Note the periodic boimdary condition 
implemented by 'coin flip' added fo fhe pafh graph shiff 
operafor fo mainfain unifarify. 


B. Interacting quantum walks 

We now pass to the analysis of the interacting quantum 
walks at the core of our invesfigafion. Given a quanfum 
walk on a simple graph G, we associate a qubit with 
each one of ifs n verfices. Then, at each step of fhe 
walk, a sequence of controlled-Z operators are performed, 
each wifh fhe coin as confrol and fhe qubif af one of fhe 
verfices as fargef. 

Schemes where the walk evolution is history 
dependent are explored in Il20l . These are similar to 
our interacting walk scheme only in that they associate 


qubits with vertices, but this is so as to record where the 
walker has been so far on fhe graph. 

Formally, fhe fofal sysfem's Hilberf space is 

Hiw = Hw<g(C2)®", (5) 

and each sfep of fhe evolufion is augmenfed wifh fhe 
applicafion of fhe Z operafor, defined as 


= \i) (Ap Z) Zc,i, ( 6 ) 

where Zc,i = [|0) (0|p (g 1* + |1) (1|(^ (g is a 

controlled-Z operation with the coin as a control and 
the T*' particle of a given laffice as fhe targef, is the 
z-Pauli operator for fhe qubif at vertex i of fhe graph, 
and is ifs idenfify operafor. The evolufion operafor of 
fhe inferaefing walk is fhen 


U = 


'n—1 


Y[t. 


{Sg g 1 g)(1p g G g 1 g)i 


(7) 


for Ip, 1(3 fhe idenfify operators on walker position 
and vertex qubit spaces respectively. In other words, 
at each step of fhe walk, fhe confrolled-Z operafor is 
applied with the coin qubit as control, and with the vertex 
qubit on the walker's current position as the target. We 
study the dynamics resulting from the application of this 
propagator with the aim of generafing enfangled sfafes 
of fhe qubifs af fhe verfices of a given graph. To fhis 
end, we consider fhe reduced sfate of fhe verfex qubif 
sysfem pcit) = Trp,c(|V^(f)) (V'(f)liw) |V’(f))iw = 
(7* |'0(O))iw sfafe of fhe overall system after the 
walker has made t steps, and investigate the (structure of) 
enfanglemenf shared by fhe qubifs af fhe verfices of fhe 
graph af hand. Clearly enough, Trp( 3 ( ) sfands for fhe 
partial trace over the walker coin and position degrees of 
freedom. 

If can be seen fhaf we have an effective control knob 
over the evolution of fhe walker in fhe form of fhe coin 
operator C, which, as previously stated, can be any 
operator in SU(2). Up to an irrelevant global phase, the 
most general form of coin operafor reads 


G(0,bi,b2) 


/g-5(71+72) cos I -e5(72-7i) sin 

V e5(7i-72) sin | e5(7i+72) cos § / ’ 

( 8 ) 


which leads fo fhe parameferisafion U U{9,4>i,4‘2) = 
nr=o^ -2i {Sg g 1 g)(1p g G(6>, <^ 1 , (^ 2 ) g 1 g)- Needless 
fo say, fhe effeef of fhe evolufion induced by fhe walk 
operafor U{d, 4>i,(j)2) depends on fhe initial sfafe of fhe 
vertex qubit system. The particular choice of coin-particle 
unitary evolution suggests the preparation for 


the vertex qubits, where |-l-) = (|0) -F |l))/-\/2 is the +1 
eigenstate of ax- In our analysis, we consider fhe walk 
initial sfafe |0) |0) S T-Lp® 'Hg> where in fhe pafh graph, 
posifion 0 corresponds fo fhe leffmosf verfex. 
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FIG. 2. Von Neumann entropy of the reduced state of the vertex qubits, the coin qubit, and the walker system (topmost to lowermost 
curve respectively). We plot Et against the discrete evolution time for a cyclic lattice of four sites [(a)-(d)] and an equal-size linear 
one [(e)-(h)]. We have used the coin parameters {9, 4>i, 412 ) £ {(37r/20, 0, 77r/20), (tt/S, 0, vr/S), (7r/4, 0, 27r/5), (27r/5, 0, Stt/IO)}. 
As the time t only takes discrete values, the continuous lines are drawn as guidance to the eye. 


III. ANALYSIS OF DYNAMICS AND ENTANGLEMENT 
SHARING STRUCTURE 

A. Entanglement across bipartitions of the system 

We have investigated entanglement in the various 
bipartitions of the system. The state of the overall system 
is pure, so we will use the von Neumann entropy to 
quantify enfanglemenf across a given biparfition. This is 
defined as 

=-Tr[pslog2Ps], (9) 

with ps the reduced density matrix of one of the 
subsystems s — P, C, G of a given bipartition. In this 
way we have considered the entanglement between 
the walker and the subsystem consishng of the coin 
and the vertex qubits. In this case then ps = 
Trc,G[|V’(f)) (V'(f)liw] G = Ln {G = Cn) for a 
linear (cyclic) graph. The results are shown in Fig. 0 
(blue (middle) curves in each panel) for laftices of both 
configurahons, n = 4 sites, and four different sets of 
coin-operation parameters. Quantitatively, the amoimt 
of entanglement is similar for the two configurations, 
and rises almost to its maximum value of 2 edits 
after a short time (the linear configuration requiring a 
relatively longer time), regardless of fhe coin-operafion 
parameters. However, a linear lattice sets a steadier 
degree of enfanglemenf in fhe time-asympfotic limif: fhe 
blue (middle) curves in Fig.|^(e)-(h) display negligible 
flucfuafions around fhe asympfofic value. 

A similar behavior is observed for the amount of 
enfanglemenf befween fhe coin's degree of freedom and 
the subsystem consisting of the walker and vertex qubits, 
i.e. for Ps = (^/’(f)liw] [cf. green (lowermosf) 


curves in Fig.|^. Obviously, in this case the maximum 
degree of enfanglemenf fhaf can be achieved is 1 ebif, as 
the coin is a qubit. However, the temporal behavior of 
enfanglemenf is much more unsfeady fhan fhaf befween 
the walker and the rest of fhe sysfem. 

Finally, when focusing on fhe entanglement between 
the vertex qubits and the joint system composed 
of fhe walker and coin, we find thaf the cyclic 
lattice configuration resembles closely the entanglement 
between the coin and the rest of the system. As the 
dimension of the vertex qubit system is 2^ = 16, the 
maximum amount of entanglement that can be shared 
across such a biparhtion is 4 edits, which is never 
achieved for the coin parameters used in our simulations. 
Our analysis demonstrates the triparhte nature of the 
walker-coin-lattice state (owing to the inseparable nature 
of the state of every biparfifion being considered ouf of 
such a grouping), as well as suggesting the negligibility 
of the entanglement between the walker and coin degrees 
of freedom in fhe cyclic configurafion. 

When assessing the entanglement brought about by 
reductions of fhe overall sysfem, we find fhaf fhe walker 
is only very weakly enfangled with the coin, in the 
case of a cyclic laftice Ill4l . On the other hand, a more 
substantial degree of enfanglemenf can be found in fhe 
linear configuration. This can be seen by considering 
the reduced state ppc = Trc^ [|'!/;(f)) ('0(f)|j^] and 
evaluating its degree of entanglement using the 
logarithmic negativity IT^ 

= max[0,log2 ||ppclli], (10) 

which is based on the negativity of parfial fransposifion 
(NPT) criterion |[T0l ITT] . Here ppQ is the partial 
transposition of ppc wifh respect to the coin subsystem. 
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FIG. 3. Logarithmic negativity of the walker-coin state ppc 
plotted against the discrete interaction time t for the four sefs 
of coin paramefers considered in Fig. (plotted in the order 
listed in Fig. with circles, squares, diamonds, and triangles 
respectively) and a cyclic lattice of n = 4 sites. We plot only the 
non-zero values of Et^. 



FIG. 4. Logarithmic negativity of fhe walker-coin stafe ppc 
plotted against the discrete interaction time t for the four sefs 
of coin parameters considered in Figs. and a linear lattice 
of n = 4 sites. We plot only the non-zero values of Et^. As 
the evolution time is discrete the continuous lines should be 
interpreted only as guidances to the eye. 


and ||X||i 


Tr 




is the Schatten 1-norm defined 


on a linear operator X. The logarithmic negativity 
was proven to be a proper entanglement monotone in 
Ref. 11311 . The results of such calculations are presented 
in Figs. 1^ and which show for the cyclic and 
linear lattice respectively, and all the coin-operation sets 
of parameters considered in our simulations [cf. Figs.|^. 


B. Entanglement between vertex qubits 


When looking into the entanglement actually shared 
by the vertex qubits, the problem of choosing a 
suitable quantifier of multipartite entanglement arises. 
While no rmique indicator is currently at hand, useful 
characterisations can be made using situation-specific 
figures of merit. One is provided by the n-partite 
concurrence, which is an extension of the widely used 
bipartite concurrence IITSi . For an n-qubit density matrix 


p this is defined as 


C{p) = max 


2 " 


0, \Ai - ^ , 

1=2 


( 11 ) 


where Ai > A 2 > • • • > A 2 " are the eigenvalues of the 
operator 

P = pSyP*Sy. ( 12 ) 

with Sy = 0^1 by, p* the conjugate of p, and ay the 
y-Pauli matrix. This measure is appealing due to its 
computationally handy form and its sensitivity to global 
entanglement. In fact, it is exactly zero if any qubit 
is separable from the rest of a system. Its non-nullity, 
though, is only a sufficient condition for multipartite 
entanglement as there are multipartite entangled states 
with vanishing n-concurrence, such as the n-qubit W- 
states cni. With such proviso, the outcomes of the 
evaluation of C(/5 g) are shown in Fig.|^ where only the 
values associated with the linear lattice configuration are 
shown, those of the cyclic one being identically null. The 
four-partite entanglement within the lattice manifests 
only at specific instants of the evolution, and remains at 
low values for most of the coin-parameter set. 

The situation can be changed rather significantly by 
exploiting the entanglement between the walker-coin 
subsystem and the vertex qubits. 

By post-selecting the state of the vertex qubits upon 
projechon of the state of the other systems, we have, 
in principle, the possibility to enhance the degree of 
multiparfite entanglement shared by the particles at 
the vertices of the lattice being considered. In order to 
provide evidence of such a possibility, here we address 
the case corresponding to tracing out the state of the 
walker and projecting the coin onto the state 

|S)p = cos p|0)p-Fe"*‘^ sinp I!)<-,. (13) 

An extensive numerical exploration of the behavior of the 
four-partite concurrence performed by taking p G [ 0 , tt] 



FIG. 5. Four-partite concurrence of the state of the vertex 
qubits for the four sets of coin-operation parameters used in 
our simulations. 
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FIG. 6. Four-partite concurrence of the conditional state of the vertex qubits obtained by discarding the state of the walker and 
projecting the coin onto the state corresponding to = 0 and p = 0. We have considered the four sets of coin-operation parameters 
used in our simulations. 


and V G [ 0 , 7 r/ 2 ], and using the four coin-parameter sets 
used in the simulations studied herein shows that, while 
such a procedure is ineffective for a cyclic lattice, it works 
well for a linear one. The maximum of C {pl^ ) is found for 
projections of the coin onto its computahonal states, i.e. 
for 1 / — 0 with p = 0 or 7 r/ 2 . The temporal behavior of the 
four-partite concurrence for such projections performed 
on the linear lattice configuration is shown in Figs. 1^ 
and 1 ^ which demonstrate the improvement achieved 
upon implementation of the conditional strategy, and 
the qualitatively similar behavior observed for the two 
projections being considered. 

The remainder of our analysis aims at shedding light 
on the actual form of the entangled states that we achieve. 


C. Trace distance from entangled states 

In order to gather informahon on the form of fhe 
state produced by the interacting quantum walk protocol 
introduced here, we quantify the closeness of pa to a 
number of "archetypal" multipartite entangled states. 
We use the quantity A(p, u) = 1 — <5(p, u), which is 
defined in terms of the trace distance between two generic 


density matrices p, a, 

= ^l|p-cr||i = (14) 

i 

where are the eigenvalues of the operator p — a. We 
consider GHZ and W states IITOl , as well as graph states 
corresponding to the lattice arrangements at hand in our 
analysis. Graph states are generalisations of cluster states, 
which are the key ingredient in measurement based 
quantum computation IT^ . Given a graph G = (G, E), 
the corresponding graph state is defined as 

|G)= n (15) 

where is the controlled-Z operator with qubit i as 
control and qubit j as target, and |-|-) = (| 0 ) -F |l))/-\/ 2 . 

We begin our characterisahon of the qubit state 
resulting from the interacting walk by calculahng its 
trace distance from GHZ, W and corresponding graph 
stafes |G„),|L„). We see from Fig. where we 
plot the quanhty A = maxg, 4 ,^,t ^ {pcit), pe) (with 
PE = |GHZ„) (GHZ„|, |W„) (W„|, and |G) (G| for G = 
G„, Ln), fhaf fhe quantum walk on the four-vertex cycle 
graph generates a perfect IG 4 ) (G 4 I state at some instant 
of time. To inveshgate this further, in Fig. we plot 
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FIG. 7. Four-partite concurrence of the conditional state of the vertex qubits obtained by discarding the state of the walker and 
projecting the coin onto the state corresponding to = 0 and ^ = tt/2. We have considered the four sets of coin-operation 
parameters used in our simulations. 


the evolution of A{pc^, IC' 4 ) (C 4 I) for the sets of coin 
parameters corresponding to the highest values achieved 
during the first 100 time steps. The highest value attained 
is exactly 1.0, and is found at f = 24 using C(|, 0, f). We 
also appreciate that at no point in the first 100 steps is the 
state of the vertex qubit system close to a W or GHZ state, 
regardless of the coin operator and graph topology. From 
this it is also clear that the linear lattice configuration 
performs markedly worse than the cyclic one. 


IV. CONCLUSION 

We have investigated a scheme for generating 
multipartite entanglement which we call the interacting 
quantum walk. As described above, this is where a 
quantum walker evolves on a graph and interacts with 
qubits located on the vertices via a phase shift interaction. 
We have performed a computer simulation of the walker 
and vertex qubit system dynamics, and have shown 
that the scheme can be used to generate multipartite 
entanglement between the qubits, as well as between 
the qubits and the walker and coin. In particular, we 
have shown that it can be used to prepare 4 qubit cyclic 
cluster states, of use in measurement based quantum 
computation. 

While here we have only considered linear and cyclic 


graph structures, a natural extension would be to 
consider different graph topologies and the effect they 
have on the amount of entanglement generated by the 
system. It is notable that the states generated by the 
interactions are never close to GHZ or W states. It would 
be interesting to find which graphs cause these states to 
be generated. 

A pre-requisite for performing measurement based 
quantum computing is a reliable method for generating 
large graph states. In flTl , Browne and Rudolph discuss 
"fusion" operators for chaining together small graph 
states to generate larger ones: an interesting direction of 
research would be to adapt these methods to situations 
where different clusters, each generated through the 
protocol presented here, are joined together via effective 
fusion gates. 

Finally, an exciting future direction would be to 
implement the interacting quantum walk in a laboratory 
setting. Experimental implementations of quantum 
walks are becoming more mature, and it could 
potentially be possible to engineer the controlled phase 
interaction required for our protocol. In parhcular, 
the use of superlattices in ultra-cold atomic settings, a 
scenario where quantum walks have been controllably 
implemented, could be a feasible avenue to explore IITHlI . 
This is particularly interesting in light of the fact 
that schemes for the generation of cluster states in 



















































7 



FIG. 8. Comparison of minimum achieved trace distance from states of interest, A for (a) pE ~ |GHZ]v) (GHZiv|, (b) pe = 
|Gjv) (GivI , \Ln) {Ln\, and (c) |Wjv) (Wiv|. 
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